arXiv:nucl-th/9701027vl 16 Jan 1997 


Three-potential formalism for the three-body Coulomb scattering 

problem * 


Z. Papp 

Institute of Nuclear Research of the Hungarian Academy of Sciences, 
P.O. Box 51, H-fOOl Debrecen, Hungary 
(8 October 1996) 


Abstract 


We propose a three-potential formalism for the three-body Coulomb scat¬ 
tering problem. The corresponding integral equations are mathematically 
well-behaved and can succesfully be solved by the Coulomb-Sturmian sepa¬ 
rable expansion method. The results show perfect agreements with existing 
low-energy n — d and p — d scattering calculations. 

PACS number(s): 21.45.+v, 03.65.Nk, 02.30.Rz, 02.60.Nm 


Typeset using REVTpX 


*to be published in Phys. Rev. C 55, March (1997) 


1 



I. INTRODUCTION 


Since the Faddeev equations are the fundamental equations of the three-body problems 
their solutions are of central interest in many branches of physics. This is especially true in 
nuclear physics because three-body calculations serve as a distinguished tool for studying 
the fundamental nucleon-nucleon interactions. A general interaction may have a local or 
non-local short-range part and a long-range Coulomb part. The solution of the Faddeev 
equations with such type of potentials is not an easy job, especially the Coulomb interaction 
and the scattering dynamics make the procedure very complicated. There exits extensive 
literature on the subject (see, e.g., Refs. jl-T] and references therein) so we restrict ourselves 
only to practical approaches. 

There are two genuinely different approaches in the practical Faddeev calculations that 
in some extent can handle Coulomb-like interactions in scattering-state problems. One of 
them is based on the solution of the configuration-space differential equations using the 
asymptotic boundary conditions ||. In the other approach, in order that the standard 
techniques could be applied, the long-range Coulomb potential is screened, and then, as 
the screened Coulomb potential goes to the unscreened one, a renormalization procedure is 
applied P|. However, in spite of the rapid development we have witnessed in the past few 
years, only limited solutions below or above the breakup threshold are available yet (see, 
e.g. Refs. @-§§. 

An another approach to the nuclear three-body problem with Coulomb interaction were 
derived along the two-potential formalism. The first, and formally exact, approach was 
proposed by Noble [fj. In this formulation all the Coulomb interactions were included in 
’’free” Green’s operator. Thus the corresponding Faddeev-Noble equations are mathemati¬ 
cally well-behaved and in the absence of Coulomb interaction they fall back to the standard 
equations. However, the associated Green’s operator is not known, so this formalism is not 
suitable for practical calculations. 

The aim of this paper is to treat the three-body Coulomb scattering problem via the 
solution of the Faddeev-Noble integral equations. In Sec. II below we shall derive a ’’three- 
potential” formalism. We will arrive at a set of Lippmann-Schwinger and Faddeev equations 
which form a mathematically well-behaved set of integral equations. In Sec. Ill below we 
shall describe how the solution can be reached. In Sec. IV below we shall compare our 
calculations with existing benchmark n — d and p — d below-breakup scattering results. It 
is found that in all cases excellent agreement is achieved. The method therefore appear as 
a promising and efficient tool for solving the Coulomb three-body scattering problem, as 
it can be adapted to more general cases and can be extended to above-breakup Coulomb 
scattering calculations. 


II. THREE-POTENTIAL FORMALISM FOR THE THREE-BODY COULOMB 

SCATTERING PROBLEM 

The Noble’s approach, which is, in fact, a two-potential formalism, requires the knowl¬ 
edge of the complete solution of the three-body Coulomb problem. Bencze has suggested 
to replace the incalculable three-body Coulomb Green’s operator by the channel-distorted 
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Coulomb Green’s operator [ 10 ]. In this section below we will follow the derivation of Ref. [10 


but instead of neglecting the intermediate range polarization potential, we will link the 
three-body Coulomb Green’s operator to the channel-distorted Coulomb Green’s operator 
via a Lippmann-Schwinger equation. Thus we will arrive at a set of Faddev-Noble and 
Lippmann-Schwinger integral equations which are mathematically well-behaved because all 
the long-range interactions are kept in the Green’s operator. 

The Hamiltonian of a three-body system with short-range plus Coulomb two-body in¬ 
teractions reads 


H — H° + v s a + v s p + id + + Vp + , (1) 

where H° is the three-body kinetic energy operator, v a denotes the interaction in subsystem 
a and the superscript s and C stand for short-range and Coulomb, respectively. We introduce 
here the usual configuration-space Jacobi coordinates and ?/ a ; £„ is the coordinate between 
the pair (/?, 7) and r) a is the coordinate between the particle a and the center of mass of 
the pair (/J, 7). Thus the potential u Q , the interaction between the pair (/?, 7), appears as 
Va (£a) ■ 

The asymptotic Hamiltonian is defined as 

H a = H° + v s a - m£, (2) 

and the asymptotic states are the eigenstates of H a 

H a \$ a ) = E\$ a ), (3) 

where {£, a Va\&a) = (Va\Xa) {£a\ 4 >a), a product of a free motion in coordinate r] a and a bound- 
state in the two-body subsystem £ Q . 

We define two asymptotic Coulomb Hamiltonians as 

= H° + v s a + v% + Vp + v° (4) 

and 

= H° + < + <£ +1£, ( 5 ) 

where u^ is an auxiliary potential in coordinate r) a , which is required to have the asymptotic 
form 

u c ~ UZjs + Z^ 

Va 

as Va — > 00. I11 fact u% is an effective Coulomb interaction between the center of mass of 
the subsystem a (with charge Zp + Z 7 ) and the third particle (with charge Z a ). 

Let us introduce the resolvent operators: 

G(z) =(z- H)-\ ( 7 ) 

GS(z) = (z-HS)-\ ( 8 ) 
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G a (z) = (z-H Q )~ 1 . 

(9) 

The operator is Noble’s channel Coulomb Green’s operator and G a is the channel dis- 

torted Coulomb Green’s operator introduced by Bencze ldj . These operators are connected 
via the following resolvent relations: 

G(z) = C%(z)+C%(,z)V‘G(z), 

( 10 ) 

GC(z) = G a (z) + G a (z)U a GC(z), 

( 11 ) 

where V a = + v ® and U a = — u 

I 11 the potential U a the Coulomb tail of Vp + is compensated by 

the asymptotic motion U a is of short-range type, so the equation ([Tl]) 
well-behaved. The scattering states 

u^- As concerning 
is mathematically 

|*S (± >) = limi e Gj(£:±i e )|t„) 

£—>0 

( 12 ) 

satisfy the Lippmann-Schwinger equations 


|$« (±) ) = |$L ±} ) + G a (E ± i0)17 Q |^ (±) ), 

(13) 

where 


|$i ±} ) = limi£G a (E ± ie) $ a ). 

£—>0 

(14) 

In conhguration-space representation the sates 4 > i ± ' ) ) are given as 


(taVal&c?) = (ValX^) (^\ K) , 

(15) 


where {rj a \Xa) are scattering functions in the Coulomb-like potential u^- 

In (jlOj) the potential V a is of short-range type and G% contains all the Coulomb interac¬ 
tions. Now, all the formulas which exist in the conventional short-range three-body theory 
can analogously be derived, only the channel Green’s operator G a has to be replaced, a la 
Noble, by G^ throughout. One can analogously perform the Faddeev decomposition and 
for the Faddeev components |of the scattering function 

= limi eG(E a ± m)|<h Q ) = |$ Q > + £ |^ ±} ) (16) 
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one arrives at the Faddeev-Noble integral equations 

ivy*} =+ c c a (E ± iojKiv-y*)+(17) 

with a cyclic permutation in a, f3, 7 . 

The S-matrix elements of scattering processes can be obtained from the resolvent of the 
total Hamiltonian by the reduction technique [|T]J 
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( 18 ) 


S/3n,am ~ lilll lim i£e 1(B »" Eam)f {$/3n\G{E am + ie) \$ am ) . 

t —>OC £—>0 

The subscript m and n denote the m-th and n-th eigenstates of the corresponding subsys¬ 
tems, respectively. If we substitute now (|T0|) into (|T8D we can get, like in [jHj, the following 
two terms: 


C _ cl 1 ! 2 ) o( 3 ) _ cl 1 ) c( 2 ) I c( 3 ) 

^(3n,am ^/3n,am ' ^ /3n,am ^/3n,am ' ^ /3n,am ' ^/3n,am 


S { Bnlm = i im limi£e 1(E ^ Eam)t ($p n \G c a (E am + U)|$ am ) 


t— HX> £^0 


s /3n,am = b™ limiee 1(E ^ Earn)t p n \G E (E am + ie)\V^G(E am + U)|<f> am ). 

We substitute again (|ll|) into (|20| ) and the first term yields again two further terms 

Sgnam = l im lim iee^n-E^t ($ | Q a (E am + ie)|$ Qm ) 

^ ’ t—>00 £—>0 


(19) 

( 20 ) 

( 21 ) 


( 22 ) 


S?L m = lim limi£e i < E 9"- E "” ), (* e „| G a (E am + i c)U a G c ’(E am + ie)|« am >. (23) 

^ ’ t —»oo £—>0 

Making use of the properties of the resolvent operators the limits can be performed and 
we arrive at the following, physically very plausible, result. The first term, , is the 

S-matrix of a two-body single channel scattering on the potential u E 

S/3n,otm = $Pot^ninS(u^). (24) 

If u E is a pure Coulomb interaction S(u E ) falls back to the S-matrix of the Rutherford 
scattering. The second term, describes a two-body multichannel scattering on the 

potential U a 

SfU = -2mS„J(Eg n - E« m )($fc > |£/“l*2, +) )- (25) 

The third term contains of the complete three-body dynamics 

■s£«™ = -2^S(E en - £U){*£- , |V 3 |'I'! 3 . (26) 

Utilizing the properties of the Faddeev components [[12] the matrix elements in (^) can be 
rewritten in a form which is better suited for numerical calculations 

(tg-y-’i'TW) = < 27 ) 

We note, that if the Coulomb interactions are absent the whole ’’three-potential” formalism 
falls back to the conventional short-range formalism. 
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III. SOLUTION OF THE THREE-BODY INTEGRAL EQUATIONS 


To solve operator equations in quantum mechanics one needs a suitable representation 
for the operators. For solving integral equations it is especially advantageous if one uses such 
a representation where the Green’s operator is simple. The free Green’s operator takes a 
very simple form in momentum representation. This is the main reason why for the solution 
of Faddeev equations, in the presence of short-range interactions, momentum-representation 
techniques perform so successfully (see for a recent review Ref. [|13|]). Since the momentum 
representation is a continuous representation, to solve the equation one needs also some kind 
of discretization. 

For the two-body Coulomb Green’s operator there exists a Hilbert-space basis in which its 
representation is very simple, it is the Coulomb-Sturmian (CS) basis. In this representation- 
space the Coulomb Green’s operator can be given by simple and well-computable analytic 

ill, 


functions 


This basis is a countable set. If we represent the interaction term on a 


finite subset of the basis it looks like a kind of separable expansion of the potential, so the 
integral equation becomes an algebraic equation. The completeness of the basis ensures the 
convergence of the method. 

In the past few years along this idea we have developed a quantum-mechanical approx¬ 
imation method for treating Coulomb-like interactions in two-body calculations. Bound- 
and resonant-state calculations were presented first [|T|], then the method was extended to 
scattering states |K| and multichannel problems |TB[. Since only the asymptotically irrele¬ 


vant short range interaction is approximated, the correct (two-body) Coulomb asymptotics 
is guaranteed. The corresponding computer codes for solving two-body bound-, resonant- 
and scattering-state problems were also published [|T7|. 

Recently the CS separable expansion approach was applied to solving the three-body 
bound-state problem in the presence of short-range plus repulsive Coulomb interactions 
18| . The homogeneous Faddeev-Noble integral equations were solved by expanding only 


the short-range part of the interaction in a separable form while treating the long-range 
part in an exact manner. The efficiency of the method was demonstrated in benchmark 
calculations of the three-body bound-state problem without and with Coulomb interactions. 
In both cases the solution showed a rapid convergence, and, whenever a comparison were 
possible to existing results in the literature, correct predictions for the binding energies and 
wave functions were achieved. The method was also applied in realistic calculations [|IIJ . 


In subsection A below we will define the basis states in two- and three-particle Hilbert 
space. In subsection B below we recapitulate some of the most important formulas of the 
two-body problem (the details are given in Refs. (T^JT^JT^] ), while in subsection C below the 
solution of the three-body Coulomb scattering problem along the CS separable expansion 
technique is presented. 


A. Basis states 


The CS functions, which are the solutions of the 
in some angular momentum state 


genic systems |20 


momentum-space as 


Sturm-Liouville problem of hydro- 
l are defined in configuration- and 
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1/2 


(28) 


(r\nl) 


{n + 2l + l)\ 


(2br) l+1 e~ br L^ +1 (2br) 


and 


(p[nl) 


2 i+3 / 2 l\(n + l + l)Vrd 

\J vr(n + 21 + 1)! 

6(26p) ;+1 i+1 / p 2 -& 2 \ 

( p 2 + 6 2)2i+2 » \j,2 + fe 2 J ’ 


(29) 


respectively, and n = 0,1,2,.... Here, L and G represent the Laguerre and Gegenbauer 
polynomials, respectively, and b relates to the energy in the Sturm-Liouville equation. We 
take b as a fixed real parameter, thus working with energy-independent bound state CS 
functions. In an angular momentum subspace they form a complete set 


N 


1 = Jim W \nl)(nl\ = Jim 1 


71=o 


N- 


N, 


(30) 


where | nl) in configuration-space representation reads 

(r\nl) = -Mnl). (31) 

r 

The three-body Hilbert space is a direct sum of two-body Hilbert spaces. Thus, the ap¬ 
propriate basis in angular momentum representation (omitting the explicit spin and isospin 
dependence from onr notation) should be defined as a the direct product 

\7ivl\) a = \nl) a <g) \v\) a , (n, v — 0,1, 2,...), (32) 

with the CS states from Eq. (|28| ) or Eq. (p9|). Here l and A denote the angular momenta 
of the two-body pair (/3, 7 ) and of the third particle a relative to the center of mass of the 
pair, respectively. Now the completeness relation takes the form (with angular momentum 
summation implicitly included) 

N 

1 = lim 'Y' \nul\) a a {nulX\ = lim 1^ (33) 

TV—>-oo z ' N —>00 

n,v=() 

where the configuration-space representation in terms of Jacobi coordinates and rj a reads: 

{^aVa\nvlX)a = (£aVa\ nvlX) a . (34) 

^aVa 

It should be noted that in the three-particle Hilbert space we can introduce three equivalent 
bases which belong to fragmentation a, [3 and 7 . 
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B. Coulomb-Sturmian separable expansion in two-body problems 


Let us study a two-potential case of short-range plus Coulomb-like interactions 

v i = vf + v c (35) 

and consider the inhomogeneous Lippmann-Schwinger equation for the scattering state \ipi) 
in some partial wave l 

\^i) = \4>i) + 9? (#KI V’i)- (36) 

Here |0p) is the regular Coulomb function, gf(E) is the two-body Coulomb Green’s operator 

sf(-E) = (E - hi - v c )-' (37) 

with the free Hamiltonian denoted by /i°. We make the following approximation on Eq. (|36|) 

\4>i) = W?) + 9 ?(E)l N vfl N \^i), (38) 


i.e. we approximate the short-range potential vf by a separable form 

N _ _ 

v s i w | ill) uf (n’l\ (39) 

n,n'=0 

where 


iL? = (nl\vf\n r l). 


(40) 


Multiplied with the CS states (nl from the left, Eq. (p8|) turns into 
equations for the wave-function coefficients 'ijj, = ( nl\ipi) 

a linear system of 


(41) 

where 


= ( nl \rf) 

(42) 

and 


9 ? n jE) = {nl\g?(E)\ril). 

(43) 


While the matrix elements of the potential may be evaluated (numerically) for any given 
short-range potential either in configuration or in momentum space, the matrix elements 
([43|) and the overlap (|4^) can be calculated analytically |T4|; the corresponding computer 
code is available from Ref. Il7l|. This fact then also allows to calculate the matrix elements 


of the full Green’s operator in the whole complex plane, 


gjz) = ((g?(z))- 1 -vtr\ (44) 

this will be needed later on in the solution of the three-body problem with charged particles. 
Of course, bound-state solutions can also be generated by solving the homogeneous version 
of Eq. ( |4l|) . 










C. Coulomb-Sturmian separable expansion approach to three-body Coulomb 

scattering problems 

In the set of Faddeev-Noble equations (|T7|) we make the following approximation: 

l« = <WO + (45) 

i.e. we approximate the short-range potential v s a in the three-body Hilbert space by a sepa¬ 
rable form 


N 

V s a « \nvix) a v s a0 0 (n'v'l'X' 

n,is,n' ,i ;'=0 


(46) 


where 


= (! - M a {nvlX\v s a \n'v'l'X') 0 . (47) 

In the ket and bra states are belonging to different fragmentations depending on the 
environments of the potential operators in the equations. 

Muliplied with the CS states a (nulX\ from the left, Eqs. (^) turn into a linear system 
of equations for the coefficients of the Faddeev components A nv — a (nz4A|^ Q; ): 


[(G 0 ) -1 - v s ]ijj = $ 


c 


(48) 


with 


££w.A'„nV = <W a {nvlX\G^\n'v'l'X') 


C I 


(49) 


and 


ib c 

—laXctTliy 


= n (ni4A|4j). 


(50) 


Notice that the matrix elements of the Green’s operator are needed only between the same 
partition a whereas the matrix elements of the potentials occur only between different 
partitions a and (3. The latter may again be evaluated numerically either in configuration 
or momentum space by making use of the transformation of Jacobi coordinates |2T|. 

Unfortunately neither the matrix elements (f|9|) nor the overlap ([50]) are known. However, 
Eqs. (HU) and (0), which are, in fact, two-body Lippmann-Schwinger equations, link them to 
relatively simpler quantities. If we perform again the separable approximation on potential 
U a with the help of the formal solution of ([TT]) we may now express the inverse matrix 
(Gf(£0 )-'as 




(51) 


where 


Qj, ot \ a nv,V a \' lx n'v' ~ a{nvlX\Ga\n'v'l'X!) c 


(52) 
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and 


= a (nvl\\U‘\nVl'\%. 


(53) 


In a similar way, with the help of the formal solution of (|T3j) we get 

iS = KSar 1 -nT 1 (s«r 1 i„ 

where 

= a{nvl\\$ a ). 


(54) 


(55) 


The state |$ a ), in fact, is a product of a two-body bound-state wave function in co¬ 
ordinate and a two-body scattering-state wave function in coordinate rj a . Their CS 
representations are known from the two-particle case of the previous section [cf. Eq. (|4l|) 1. 

For the calculation of the matrix elements in Eq. ( |5| ) we proceed in a similar way as in the 
case of three-body bound-states ||T8|. Since in H a of Eq. (||) we can write the three-particle 
free Hamiltonians as a sum of two-particle free Hamiltonians 


= hi + hi, 


(56) 


the Hamiltonian H a appears as a sum of two Hamiltonians acting on different coordinates 

= h$ a + h Va , (57) 

with li£ a = + v a(£a) + v a'(£a) and h Va = h® a + u^(rj a ), which, of course, commute. Thus 


we can apply the convolution theorem [B2 


G a (z) = {z- h ia - h Va ) 1 = £ dw(z -w- h v J 

x (w ~ %J _1 . 


-i 


(58) 


Here the contour C should encircle, in positive direction, the spectrum of h^ a without pen¬ 
etrating into the spectrum of h Va . For scattering-state energies at real z these singularities 
overlap. To find the correct path one should take the z = E + ie case with hnite e. Now 
the condition on C can easily be fulhlled and then, one should take the £ —> 0 limit allowing 
only analytic deformation for the contour C [see Fig. [l| . 

After sandwiching the above Green’s operator between the CS states, the integral in Eq. 
(j58[) appears in the form 

G-l a X ol nv,l' a X' a n'v'( E + iO) 

1 


(59) 


= f c dw + iO~w- h Va ) 

x a (nl\(w - ^ a ) _1 |n'/') a , ; 


where both matrix elements occurring in the integrand are known from the two-particle case 
[cf. Eq. (||)]. 
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IV. TESTS OF THE METHOD 


In this section we demonstrate the performance of the method in calculations of three- 
body short-range and Coulombic scattering phase shifts at energies below the breakup 
threshold. We have selected cases that serve as benchmarks for various three-body scat¬ 
tering calculations. As an example we take a model three-nucleon problem with s-wave 
Malflet-Tjon (MT) I-III potential, acting in singlet and triplet states, as parametrized in 
Ref. H . We have calculated quartet and doublet n — cl and p — d phase shifts and compare 
them to the results of the configuration-space Faddeev calculations of Ref. || . 

Before presenting the final results, let us demonstrate the convergence of the results for 
scattering phase shifts at various energies. We take extreme cases, one is at very low energy, 
an other one is jut below the breakup threshold, the third one is in between. We select 
two-channel doublet n — d and p — d cases, because this case is more complicated than 
the one-channel quartet case. Tables | and [TTj show that convergence up to 4 significant 
digits can comfortably achieved with N = 30 terms applied for n and v in the separable 
expansion. Remarkably, the speed of convergence is everywhere similar, irrespective of 
energy and whether or not Coulomb forces are present. 

In Tables |Hj and [TV] we compare our converged results to the configuration-space Faddeev 
calculations of the Los Alamos-Iowa group [§j. We can report perfect agreements in all cases. 

In Tables [TJ and E3 and also in Ref. [6|] the Coulomb modified phase shift from short- 
range plus polarization potential (6 c,ps ) were presented, and the corresponding values were 
used in the calculation of the scattering lengths a p f s . However, since our three-potential 
formalism allows a unique separation of these two effects we can also calculate the Coulomb 
plus polarization modified short-range phase shift 5 cp,s and the corresponding scattering 
lengths a c p /. Theoretically the scattering length should be calculated from 5 cp,s since a p % s 
is minus infinite [^3|. In practical calculations the extrapolation to zero energy were made 
from higher energy and this minus infinity limit were not seen. Careful analyses of low- 
energy p — d calculations indicated that in these condition a p % s is a good approximation to 
a pd S 0IHH. We have calculated the scattering length from both phase shifts using the 


formula of Ref. 23 


a = lim a(k) = lim 

k^o k^o kC 2 (k ) 


(60) 


where k is the wavenumber and 


C 2 (k ) = 


2,717/ 


e 2nr, _ 


( 61 ) 


with r/ = me 2 /hk, m beeing the reduced mass. From a{k ) values correspondig to energies 
down to 0.01 MeV in the doublet case and to 0.001 MeV in the quartet case, where the phase 
shift still can reliably be calculated, we extrapolated to zero energy. The results are given 
in Table 0. We can see that the effect of the polarization potential can realy be neglected 
in the calculation of the p — d scattering lengths. 

Besides the number of terms in the expansion there is only one parameter in the method, 
the b parameter of the basis. This should be chosen according to the range of the potential. 
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We have found that the converged results do not depend on b and for a wide range of 
reasonable values even the speed of the convergence is rather insensitive to the choice of b. 
In all calculations presented the same value was applied (6 = lfm -1 ). 


V. CONCLUSION 


We have suggested a three-potential formalism for treating the three-body Coulomb 
scattering problem. In absence of Coulomb interactions the formalism falls back to the usual 
short-range formalism. According to the three-potential picture the three-body Coulomb 
scattering starts with a two-body single channel Coulomb scattering, then it goes over to a 
two-body multichannel scattering on the intermediate-range polarization potential, finally 
comes the three-body scattering due to the short-range potentials. This formalism preserves 
the mathematical correctness of Noble’s approach, and along the idea of channel distorted 
formalism of Bencze, without neglecting important terms, gives solvable equations. 

These ’’solvable” equations are certainly too complicated for most of the numerical meth¬ 
ods available in the literature, but the Coulomb-Sturmian separable expansion method can 
cope with them. It solves the three-body integral equations by expanding only the short- 
range part of the interaction in a separable form on a Hilbert-space basis while treating 
the long-range part in an exact manner via a proper integral representation of the three- 
body channel distorted Coulomb Green’s operator. As a consequence the method has good 
convergence properties and can in practice be made arbitrarily accurate by employing an 
increasing number of terms in the expansion. The usage of the Coulomb-Sturmian basis 
is essential as it allows an exact analytic representation of the two-body Green’s operator, 
and thus the contour integral for the channel distorted Coulomb Green’s operator can be 
calculated also in the practice. 

We have presented below-breakup calculations and got perfect agreements with existing 
benchmark results. We have observed a fast convergence with respect to increasing the 
number of terms in the expansion. Using high rank expansion we got very accurate results. 
The quality of the convergence is practically the same in short-range and in Coulomb case, 
Coulombic calculations need only roughly 30% more computer time. Although the example 
presented here is rather simple, but not unrealistic, the method can handle more complicated 
potentials, as were demonstrated in bound state calculations [ |T8| , |T9[ | . 

Certainly, the toughest problems in nonrclativistic three-body scattering are the above- 
breakup calculations with Coulomb interactions. In this respect, the method presented here, 
is very promising since the equations used are mathematically well-behaved also for this case. 
The extension of the contour integral for above-breakup energies is straightforward. The only 
foreseeable problem is that the interaction volume is much bigger and one needs much higher 
terms in the expansion. Indeed, test calculations show that for energies just a little bit above 
the breakup threshold with terms up to IV = 34 we can reach acceptable convergence and 
good agreements with benchmark results ^5f, but the method fails to reach convergence 


for higher energies. This indicates that the mathematical formulation is correct, but the 
available computing power is not sufficient. So, the method needs some more polishing and 
we have to think a little bit further. 
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TABLES 


TABLE I. Convergence of the 2 fi n d phase shifts for three-nucleon system interacting via the 
MT I-III potential at various energies, with increasing basis for the separable expansion. N denotes 
the maximum number of basis states employed for n and u. The phase shifts are in degrees. 


N 

0.1 MeV 

2 find 

1.0 MeV 

2.18 MeV 

10 

-4.0908 

-20.704 

-32.896 

11 

-3.6803 

-20.330 

-33.596 

12 

-3.8876 

-20.450 

-33.599 

13 

-3.4583 

-20.630 

-33.901 

14 

-3.5189 

-20.562 

-33.963 

15 

-3.3917 

-20.629 

-33.834 

16 

-3.3434 

-20.653 

-33.723 

17 

-3.3293 

-20.636 

-33.594 

18 

-3.2914 

-20.660 

-33.525 

19 

-3.2940 

-20.652 

-33.341 

20 

-3.2764 

-20.654 

-33.282 

21 

-3.2783 

-20.656 

-33.267 

22 

-3.2709 

-20.653 

-33.290 

23 

-3.2715 

-20.655 

-33.326 

24 

-3.2688 

-20.654 

-33.377 

25 

-3.2688 

-20.654 

-33.424 

26 

-3.2681 

-20.654 

-33.460 

27 

-3.2678 

-20.654 

-33.483 

28 

-3.2678 

-20.654 

-33.491 

29 

-3.2676 

-20.654 

-33.486 

30 

-3.2676 

-20.654 

-33.473 

31 

-3.2675 

-20.654 

-33.456 

32 

-3.2675 

-20.654 

-33.438 

33 

-3.2675 

-20.654 

-33.423 

34 

-3.2675 

-20.654 

-33.413 
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TABLE II. Convergence of the 2 5 p d phase shifts for three-nucleon system interacting via the 
MT I-III potential at various energies, with increasing basis for the separable expansion. N denotes 
the maximum number of basis states employed for n and u. The phase shifts are in degrees. 


N 

0.1 MeV 

2 fi p d 

1.0 MeV 

2.0 MeV 

10 

-0.9485 

-16.376 

-28.405 

11 

-0.7947 

-15.924 

-28.659 

12 

-0.8794 

-16.089 

-28.564 

13 

-0.6721 

-16.195 

-28.847 

14 

-0.7015 

-16.138 

-28.874 

15 

-0.6241 

-16.210 

-28.840 

16 

-0.5962 

-16.226 

-28.865 

17 

-0.5817 

-16.216 

-28.840 

18 

-0.5575 

-16.241 

-28.802 

19 

-0.5572 

-16.232 

-28.800 

20 

-0.5446 

-16.238 

-28.774 

21 

-0.5455 

-16.238 

-28.767 

22 

-0.5390 

-16.236 

-28.764 

23 

-0.5394 

-16.238 

-28.761 

24 

-0.5362 

-16.236 

-28.765 

25 

-0.5363 

-16.237 

-28.768 

26 

-0.5350 

-16.237 

-28.772 

27 

-0.5349 

-16.237 

-28.776 

28 

-0.5345 

-16.237 

-28.778 

29 

-0.5343 

-16.237 

-28.780 

30 

-0.5342 

-16.237 

-28.780 

31 

-0.5341 

-16.237 

-28.780 

32 

-0.5341 

-16.237 

-28.779 

33 

-0.5340 

-16.237 

-28.778 

34 

-0.5340 

-16.237 

-28.778 
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TABLE III. ^5 nd and 2 5 nd phase shifts for three-nucleon system interacting via the MT I-III 


potential at various energies. 

The phase shifts are in degrees. 



E 

Ref. § 

4 5 nd 

This work 

2 

Ref. [| 

$nd 

This work 

0.001 

-2.09 

-2.092 

-0.230 

-0.229 

0.05 

-14.6 

-14.62 

-1.99 

-1.988 

0.1 

-20.4 

-20.44 

-3.28 

-3.267 

0.2 

-28.3 

-28.29 

-5.68 

-5.670 

0.3 

-34.0 

-33.98 

-7.95 

-7.944 

0.4 

-38.5 

-38.54 

-10.1 

-10.09 

0.5 

-42.4 

-42.37 

-12.1 

-12.12 

0.6 

-45.7 

-45.69 

-14.0 

-14.03 

0.7 

-48.6 

-48.63 

-15.8 

-15.83 

0.8 

-51.2 

-51.27 

-17.5 

-17.53 

0.9 

-53.6 

-53.66 

-19.1 

-19.13 

1.0 

-55.8 

-55.86 

-20.7 

-20.65 

1.633 

-66.7 

-66.72 

-28.6 

-28.60 

2.180 

-73.6 

-73.6 

-33.6 

-33.4 


TABLE IV. 4 S pd and 2 6 pd 
potential at various energies. 

phase shifts for three-nucleon system interacting via 
The phase shifts are in degrees. 

the MT I-III 

E 

Ref. § 

%d 

This work 

2 s pd 

Ref. [| 

This work 

0.001 

0.0 

0.0 

0.0 

0.0 

0.05 

-2.69 

-2.694 

-0.113 

-0.112 

0.1 

-7.46 

-7.458 

-0.537 

-0.534 

0.2 

-15.6 

-15.56 

-1.96 

-1.949 

0.3 

-21.9 

-21.86 

-3.73 

-3.720 

0.4 

-27.0 

-27.00 

-5.62 

-5.612 

0.5 

-31.3 

-31.34 

-7.53 

-7.520 

0.6 

-35.1 

-35.11 

-9.40 

-9.394 

0.7 

-38.4 

-38.43 

-11.2 

-11.21 

0.8 

-41.4 

-41.40 

-13.0 

-12.96 

0.9 

-44.1 

-44.09 

-14.6 

-14.63 

1.0 

-46.5 

-46.55 

-16.2 

-16.24 

0.667 

-37.3 

-37.37 

-10.6 

-10.62 

1.333 

-53.5 

-53.49 

-21.1 

-21.08 

2.0 

-63.8 

-63.74 

-28.8 

-28.78 
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TABLE V. p 
potential. 

— d scattering lengths for three-nucleon system interacting 

via the MT I-III 


This work 

Ref. |H 

4 rC,pS 

° P d 

13.76 

13.8 

4 rcp,s 
°pd 

13.79 


2 :C,ps 
°pd 

0.161 

0.17 

2 rcp,s 
°pd 

0.195 
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FIGURES 



the topology of the contour should be kept. 
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